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i The fluctuation induced diamagnetism of underdoped high temperature superconductors is stud- 

ied in the framework of the Lawrence-Doniach model. By taking into account the fluctuations of 
' the phase of the order parameter only, the latter reduces to a layered XY-model describing a liquid 

t of vortices which can be either thermally excited or induced by the external magnetic field. The 

diamagnetic response is given by a current-current correlation function which is evaluated using the 
Coulomb gas analogy. Our results are then applied to recent measurements of fluctuation diamag- 
netism in underdoped YBCO. They allow to understand both the observed anomalous temperature 
dependence of the zero-field susceptibility and the two distinct regimes appearing in the magnetic 
(— i | field dependence of the magnetization. 

Q . PACS: 74.20. De, 74.25. Ha, 74.40. +k 
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I. INTRODUCTION 

S Owing to their short coherence length, high temperature superconductors show marked deviations from the mean- 
field behaviour that describes rather well the behaviour of conventional superconductors. The underdoped regime 
j of the various cuprates is particularly interesting, given their pronounced anisotropy and the low density of charge 
c~j . carriers. In this region of the phase diagram, fluctuations are strongly enhanced and are manifest already well above 
the critical temperature T c . This fluctuation regime gives rise to various unusual phenomena,!!!! such as an anomalous 
temperature dependence of the Knight shift, NMR relaxation rate and electrical conductivity, an anomalous frequency 
dependence of infrared conductivity, as well as the formation of a pseudogap in the electronic density of states—Close 
t— ( | to T c it also allows to see true critical behaviour in quantities such as the specific heata or thermal expansion.tj 
^ ■ In the present work, we specifically address the temperature and field dependence of the diamagnetic susceptibility 
X of strongly anisotropic high temperature superconductors above T c . The influence of fluctuations on \ nas been 
studied long ago in the framework of a LandaujXJinzburg model, taking into account Gaussian fluctuations of the 
order parameter above the critical temperature.u More refined calculations based on the Lawrence-Doniach model 
have taken into account the lattice structured Recent measurements^ have shown that the Gaussian approximation 
can indeed well describe the diamagnetic fluctuations in optimally doped YBCO, whereas underdoped specimens of 
the same compound, such as YBa2Cu30g.67j show dramatic deviations from this behaviour. The fluctuation region 
where the zero-field orbital susceptibility shows appreciable values extends over a much larger temperature range than 
in the optimally doped system. Moreover, the field dependence of the magnetization shows a much more pronounced 
crossover between low and high fields than what would have been expected from Gaussian fluctuations. 
We base our calculations on an anisotropic Lawrence-Doniach (LD) functional, involving the superconducting order 
parameter field A in the presence of a vector potential that describes a homogeneous magnetic field perpendicular 
to the lattice planes. Rather than considering Gaussian fluctuations, we assume a "precursor regime" in which the 
£j , amplitude of A has already acquired a non-zero average value whereas its phase is subject to strong fluctuations 
inhibiting long range superconducting order. In this context the LD- functional reduces to an anisotropic layered XY- 
model. The relevant thermal excitations of such a system are the phase field singularities which manifest themselves 
as vortices and antivortices in 2D and as vortex loops in 3D. The applied magnetic field also acts on the phases by 
inducing vortex lines crossing the sample from one end to the other. The diamagnetic susceptibility, expressed in 
the usual way by a current-current correlation function, is then related to the positional correlation function of the 
vortex line elements, the static structure factor S'(q). We model S'(q) in a simple way by using various known results 
obtained either by analytic considerations based on the Coulomb gas analogy or by Monte Carlo (MC) simulations 
of the anisotropic 3D XY-model. Then we arrive at explicit expressions for the zero-field susceptibility x(T) and the 
temperature and field dependent magnetization M(T, B). They contain several material dependent parameters which 
are estimated by comparing with the experimental data on underdoped YBaCuO from Ref.6. 

This procedure unravels three main features. First, the density of the vortex line elements contributing to the 
diamagnetic response (those which are oriented in z-direction, i.e., parallel to the applied field) is thermally activated 
with a value of the activation energy that is compatible with what is found in the above mentioned MC simulations. 
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Second, in the temperature range covered by the experiments, the value of the anisotropy shows that the positions 
of the vortex line elements fluctuate strongly from one layer to the other. This points to a rather weak effective 
coupling between layers, which is compatible with the fact of being above the vortex melting, respectively the "vortex 
decoupling" line. Finally the rather sharp crossover of the magnetization M (T, B) as a function of B points to a 
subtle interplay between thermally excited vortex loops and field induced vortex lines. 

These observations allow for the following conclusions: (i) The experiments in Ref.6 can be much better accounted for 
by phase excitations than by Gaussian fluctuations of the pairing field. This is particularly manifest in the activated 
T-dependence of the susceptibility, but also in the existence of two field regimes with quite different behaviours of the 
magnetization. In this context we have, however, to admit that the very sharp crossover between these two regimes 
observed in the experiment may also be due to sample inhomogeneities as it was suggested in the experimental papers.El 
(ii) Except for the data taken close to the zero-field critical temperature, the diamagnetic response of the underdoped 
compound presented in Ref.6 seems to be "precritical" , in the sense that the relevant lengths in the lattice planes 
as well as in the perpendicular direction do not show any true critical (i.e., singular) behaviour (which would be 
supposed to belong to the 3D XY universality classa). 

In section II we develop the theoretical formalism that allows to express %(T) and M(T, B) in terms of the vortex line 
structure factor, and in section III we compare our theoretical results with the data presented in Ref.6 for underdoped 
YBCO, thereby extracting the free parameters from the experimental curves. A summary is presented in section IV. 



II. DIAMAGNETIC RESPONSE 



We discuss the orbital magnetic response of an underdoped superconductor in the London approximation to the 
Lawrence-Doniach (LD) model, i.e., in the framework of an 3D anisotropic layered XY-model in which the phase 9 
of the superconducting order parameter A is coupled to the vector potential A|| describing a homogeneous magnetic 
field B perpendicular to the lattice planes (we restrict ourselves to temperatures T > T c where the Meissner effect is 
absent, identifying thus the external and the effective internal vector potential): 

W ^ = 2^d£ / ^U||a 2 [V||^-|^A||] 2 + Jx[l-co S (^-^ +1 )]}. (1) 

n 

Here a is the lattice constant in the planes whereas d denotes the distance between two layers. Jm and J_l are the 
respective phase couplings. Their ratio 

7 2 = J > I (2) 

determines the anisotropy of the system. The XY-Hamiltonian can be obtained starting from a LD functional for the 
complex superconducting-pairing field A by keeping the amplitude of the latter constant (London approximation). 
This is a current strategyQ based on the assumption that the various precursor phenomena, observed in underdoped 
cuprates between T c and some higher temperature T* and mentioned in the introduction, are essentially due to 
fluctuations of the phase of A whereas its amplitude |A| maintains a finite mean value. Moreover, these materials 
exhibit a layered structure which is specifically taken into account by the Lawrence-Doniach approach. Expression 
(Q) is a partial continuum version of the discrete 3D XY-model, the Josephson coupling 1 — cos(#i — 9j) between 
neighbouring lattice sites in a given layer having been replaced by the phase gradient. 

The orbital magnetic response A in a finite external field is obtained by adding a small perturbation to the applied 
vector potential 

A|| -f A|| + <5A|| (3) 

and by calculating the second derivative of the free energy with respect to the perturbing field SAu. The magnetic 
susceptibility x(T, B) is then given by 

X = l Un o^~ (4) 



with 



d v< J>o k B T 
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The second term of (||) is the diamagnetic response, whereas the first term involves the current-current correlation 
function 

= Z2 £ / d2pd2p ' eiq,(r " r,) 0*(r)j*(r')> (6) 
j x (r) = VA(p) - ^A {lx (r), (7) 

Here the coordinate r means (p, nd) and the sample volume is Vl — L 2 Nd where N is the number of layers. In order 
to have a finite x, the limit C'(q — > 0) has to cancel the diamagnetic term. It will be shown below that this is indeed 
the case in our approach. We have chosen the gauge in which 

Ay (r) = 0,0). (8) 

Thus the relevant wave vector q has only a y-component denoted simply by q. 

In order to study the thermodynamic properties of the 3D anisotropic XY-model, we first recall that the Berenzinskii- 
Kosterlitz-Thouless transition occurring in the strictly 2D case is best described in terms of vortex and antivortex 
excitations. Although the 3D XY-system shows a "normal" second order transition, even when it is anisotropic, it has 
been shown that it is also possible to understand the critical behaviour of such a system in terms of vortex excitations, 
which - for topological reasons - now have to form closed loops or continuous lines crossing the whole system.Ei The 
loops are the 3D extension of the planar vortex and antivortex structure whereas the lines arise from the presence of 
an external magnetic flux penetrating into the sample in the same way as in a type II superconductor below T c . In 
the following we will use this vortex picture of the 3D XY-model in order to calculate the phase correlation function 
C(q) that determines the diamagnetic response according to Eq.(^). 

Since vortex lines are either closed (forming a loop) or extend continuously through the whole sample, we can 
characterize their structure by labelling each line by an index s and by giving its position R(s, n) in a given layer n 
which corresponds to the centre of the corresponding "pancake vortex" (see Fig.yJ). The x-component of the phase 
gradient created by all the vortex lines is then given by the same expression used in magnetostatics in order to 
calculate the magnetic field of a system of current loops and lines: 

y x 9 n (p) = d e xa p I dp 1 - . (9) 

Here e a ^ 7 is the fully antisymmetric tensor of rank 3 and the vector field K(s,r') is given by the line clement 
"tangential" to the vortex line number s at point r' = (p',n'd): 

K(s, r') = t(s, n') 5(p' - R(s, n')) z, (10) 

z being the unit vector in the z-direction and the sum in (^|) thus runs over all vortex line elements s and layers n. In 
order to make connection with the (more simple!) 2D case, it is useful to attribute a topological number t(s, n) = ±1 
to each vertical vortex line element. Its sign is chosen such that the product t(s, n) z gives the oriented "tangential" 
vector of the vortex line at that point [i.e., t = +1 (—1), when the line moves upward (downward) with respect to the 
lattice plane n]. The current correlator (|6|) is then given by two contributions: 

C(q , = I ^ |s( „-(^)X ol . (id 

The first term, stemming from the phase gradient in the current density (^) , represents the structure factor of the 
vortex line elements oriented in z-direction, given by 

5(q)= t(s,n)t(s',7i')(e lq - [R(s '" ) - R(s ''™' )1 ). (12) 

The other term in (11) comes from the second (diamagnetic) contribution to the current (0) and represents the total 
flux going through the system due to the applied field (cross correlations between V x 9 and An x are supposed to vanish 
due to the disordered structure of the vortex system above the melting temperature) . In the Appendix we show that 
the singular zero wave vector value of the correlation function S(q) compensates the second contribution to (|lTJ) so 
that the limit C(q — + 0) is well behaved. The main problem left is thus to find a suitable form for the regular part of 
the correlation function S(q). 
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First, we recall that we are interested in the fluctuation induced diamaffnetism above the zero-field transition tem- 
perature T c , a range in which recent experiments have been carried outfl For a finite field one has thus to deal with 
a region in the T-B-phase diagram laying beyond the different possible tr.aiisi|tion lines which are discussed in the 
literature, such as the vortex melting line and the vortex decoupling line.EjO Our structure faster <S(q) thus has 
to describe a disordered vortex-liquid system in which strongly "wrinkled" vortex lines and loopsBO go through the 
sample. We use the following reasoning for obtaining an approximate form of the regular contribution to S'(q): 

(i) In the extreme limit of totally decoupled layers, each plane would have to be described by a 2D XY- model 
above the critical temperature, which is frustrated when an external magnetic field is applied. Its vortex structure 
would be the one of a neutral Coulomb gas, more precisely of a mixture of a neutral two-component Coulomb gas 
(given by an equal number of thermal vortices and antivortices) and of a one-component gas (the field induced vor- 
tices) in a neutralizing background (given by the external flux). In the purely 2D case, the non-trivial part of Eq.([T^) 
thus reduces to 

fiMq) =N^t{s, 0) t(s', 0) ( e ii'[R(s.o)-R(s',o)]) 

(13) 

= NL 2 n v S c {<i). 

Here we represented the positional structure factor involving a sum over all_jzprtex and antivortex positions by the 
Coulomb gas structure factor S'c(q). We use the following approximate formEj 

Sc(q) = 2 ^ 9 ^ 2 7T T (14) 
q 2 + 2nn v q v /k B T 

wh£re ny is the areal vortex density and the "charge" qv of each vortex is related to the in-plane phase coupling Jm 

i.yH 

£#. = 27rJ||. (15) 

The expression ([l4]) yields the correct limiting behaviour of Sc (q) for q — > and should be valid for temperatures not 
too close to T c . Inserting expressioixi|l^) into Eq.([ll]) gives a similar form of the current-current correlation function 
C(q) as used by Kwon and DorseyM 

(ii) For a strong anisotropy and for the considered temperature range the effective interlayer coupling will be very 
weak. For evaluating expression (|l2"|), we split the double sum over n and n' into two parts. In the first part, we take 
n equal to n' which yields a 2D problem analogous to the one treated above. The corresponding contribution to jl^ ) 
is therefore similar to the structure factor of point vortices in a single XY-plane: 

Si(q) - NL 2 n v S c {q) (16) 

where ny is now the areal density of vertical vortex line elements. The second part 5 2 (q) of the double sum in 
( |l2| ) contains n and n' involving different layers. Here we use the fact that in the considered temperature range, the 
vortex loops and also the field induced lines have a verv-irregular shape. Indeed, the region in the T-S-plane we are 
interested in lies above any lattice layer decoupling linc£3 This fact is formulated in Ref.7, by indicating a T-domain 

IT - T r 



> (17) 



with an XY-exponent v ss 0.6, for which the lattice planes are practically decoupled. Therefore only those ones which 
are close enough to another (in a sense to be specified below) are correlated and will contribute to S^q). We first 
expand formally in expression (|l2|) the position R(s', n') of a the vortex line s' in layer n' with respect to its value in 
layer n. This gives 



S 2 (q)= t{s,n)t{s\n / )(e t * [R(s ' n) - R{s '' n)] e- t * u ( s '' n '- n '>). (18) 



Then we split the average bracket of (J18|) and factor out from the sum over n' the first exponential which pertains to 
a given layer n as well as the topological numbers t(s\n') which are equal to t(s',n) for neighbouring layers n and 
nf. This yields the previous result (|l^). We are then left with correlations between the positions R of a single vortex 
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line s' in layer n' and in layer n ^ n' . Assuming that such correlations are the same for all vortex lines s' and extend 
only over a distance £3 = n 3 d, we obtain 



"3 



5 2 (q) = 2 5i(q) V (e-*»-W) =25 1 (q)A-(q). (19) 



=1 



Here u(n) is the deviation of a given line or loop from a straight line along the z-direction. The factor 2 comes from 
the fact that the sums over n < n' and n > n' have been reduced to one such sum in (|l^) . Assuming that the vortex 

lines behave like harmonic strings with an effective stiffness given approximately by A ~ ^ = 3^2 > one finds 

"3 -. 
X(q)«^exp(-- g 2 (u 2 (n))) 

n=1 (20) 
_ 1 -exp(-k B Tn 3 dq 2 /4:X) 

~ exp(-fc B Tdq 2 /4A) - 1 

Here we used the relation (u 2 (n)) = | n d applying for harmonic deformations and yielding a geometric series. 
The final result for the regular part of C(q) in (O) then takes the form 



C(q) = ^-nW? c (q)[l + 2X(q)]. (21) 
T 

In order to evaluate x according to (|J), we have to expand C in powers of q. First, the above mentioned cancellation of 
C(q — 0) in (||) with the (first) diamagnetic contribution is fulfilled, provided that the effective charge of the Coulomb 
gas structure factor is chosen to be 

q 2 v = 2tt J,, (1 + 2n 3 ) (22) 

which is a reasonable generalization of the purely 2D result (|l5| ) , taking into account the fact that the "charges" are 
now vortex lines elements correlated over a distance ~ n 3 d along the z-direction. The cancellation of the q = term 
in Eq. ^ guarantees that the phase system of the superconductor has no stiffness above the critical temperature by 
making the limit (0) finite. Evaluating the latter leads directly to the final result for the bulk susceptibility \ (P er 
unit volume) which reads 

' / "'' 7 ' 1 +(7r 7 ) 2 n 3 (l + n 3 )d 2 ]. (23) 



(l + 2n 3 )d % L n v 

Expression ( p3[ ) will be used in the following section for interpreting the experimental data (low field susceptibility 
and field dependence of the magnetization) obtained on underdoped YBCO.H 



III. ANALYSIS OF EXPERIMENTAL MAGNETIZATION AND SUSCEPTIBILITY DATA 

In order to apply our theoretical results to underdoped YBCO, we use the following (approximate) values for the 
lattice parameters : a — 4A and d — 12A. Moreover, the molar quantities like those reported in Ref.6 are obtained by 
multiplying expression (|2^) by the molar volume Qjvf = MAa 2 d = 115 cm 3 and by a reduction factor A which accounts 
for the fact that in the materials we are interested in only a fraction of the unit cell volume actually carries the current 
densities and thus contributes to the diamagnetic response. 



A. Low-field susceptibility 

The susceptibility data of Ref.6, for a very low applied field of 0.02T, cover the range from T = 63K, which is 
just above the zero-field transition temperature, to T = 110K. Above T = 80K the diamagnetic susceptibility x is 
essentially equal to zero, a background consisting of spin susceptibility and free-electron orbital diamagnetism having 
been subtracted. For zero applied field, our expression ( p3| ) for x contains several yet undetermined parameters: 
the areal density ny of vertical vortex line elements that are thermally excited, the correlation length n 3 d and the 
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anisotropy 7. We mention that the subsequent analysis of the temperature and field dependent magnetization will 
show that the second term of (|23| ) is irrelevant at low fields. Therefore the value of 7 is unimportant and the dominant 
contribution to x(T, B — > 0) is given by in-plane correlations of the thermal vortex loop elements. Fig. 0, in which the 



data of Ref.6 are reproduced with our expression (23), shows that the observed zero-field diamagnetic susceptibility 
is almost perfectly fitted by 

X (T,B = ) = Ccxp(^) (24) 
k B T 

with E /kBT c w 22. This means that the quantity ny must be temperature dependent and obey a thermally activated 
behaviour 

n v = n ex P(~J^f) ( 25 ) 

with the factor no ~ 10 /a 2 . The number of vortex excitations in the XY- model is indeed known to show axL.ac±ivated 
T-dependence. In 2D this has recently been confirmedllJ and various authors find the same ces ult in 3IXfra using 
numerical studies. In the first case the activation energy is roughly given by Eq ~ 10 fesTcO Simulations for the 
anisoti'opic 3D case yield values that are somewhat largerEj and that depend on the structure of the corresponding 
loopJi3c3 getting larger the more planes are crossed. In a strongly anisotropic 3D XY-systpmy-jmast, of the thermal 
loops existing up to the critical temperature consist of elements parallel to the lattice planes .□I23€3£JI23 Loop segments 
perpendicular to the planes only begin to be formed above T c . Thus the activated form attributed to the vortex line 
density ny in ( pq ) can be seen as a measure for the rapidly increasing number of loops crossing two planes or even 
more. Such loops contain at least two line segments perpendicular to the planes, the rest being between two planes. 
This 3D structure may explain the relatively large value of the activation energy Eg found above. The latter is also 
compatible with an explicit expression for the loop self-energy proposed in Ref.7 assuming that the total length of the 
loop segments oriented in z-direction is on the order of two interplanar distances, the rest of the loop being oriented 
parallel to the lattice planes. The value of Eq further suggests that the quantity n 3 d measuring the extension of 
the correlations of the vortex structure along the z-direction must be of the order of a few lattice distance d (we 
will take 713 = 2 in the following). This provides an a posteriori justification for the assumptions we made when 
deriving Eq.(pl|). Finally we note that the fact that the diamagnetic response is associated with loops containing line 
elements parallel to the z-direction directly follows from the magnetostatics: the vertical loops elements are generated 
by currents flowing perpendicular to the external field B and thus responding the most sensitively to the latter. 
We can now estimate roughly the value of the density of the vortices contributing to the diamagnetic response in the 
temperature range of interest. Using the above values of Eq and no and taking into account a reduction factor of 0.15 
(see below), wc get ny ~ 10 2 /im~ 2 for T = 65K. This relatively small value will be discussed in the next section but 
we should not forget the fact that it concerns only the part of the total number of vortices which contains elements 

perpendicular to the layers. We also remark that ytiy 1 3> a, which fully justifies the continuous approach in the 
lattice planes used in the LD-action (|l]). 

B. Field dependence of the magnetization 

Here we have to deal with coexisting thermal loops and field induced vortex lines. 
The magnetic field dependence in expression ( p3| ) for \ is hidden in the density ny and, possibly, in 77,3 d, the interplanar 
correlation length. The most simple approach to deal with this situation consists in splitting ny into a thermal and 
a field induced part as follows: 

n v = nf + n{, = nf[l + z{T,B)]. (26) 
This gives naturally rise to the dimensionless variable 



/ 

^ - (27) 



n J v _ B 

^th „~>th 



$0 n\ 



Ly -1-U'l-y 



quantifying the relative importance of the two types of vortex elements. Assuming that 713 does not depend on B, 
one can integrate x(^; B) in order to obtain the magnetization per unit volume: 



M(T.B) - * ^[$ t , l„j.(l :{T.D)) + U-)-n,(l + n ^ B . r2s> 
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The corresponding molar quantity is multiplied by the additional prefactor Xfl m already mentioned at the beginning 
of this section. M depends on temperature through n'y and, possibly, through n.3. The first contribution to (|2§|), 
given by log(l + z), has two different limiting behaviours. For low fields log(l + z) <~ z and M is essentially given by 
x{T, B = 0) B, whereas for high fields the first part of M varies like log(z). The crossover between the two regimes 
takes place for z ~ 1 which describes the situation where the quantities of thermally excited and field-induced vortex 
line elements are similar to one another. Thus the corresponding field value depends crucially on the thermal vortex 
element density riy 1 showing up in the dimensionless variable z defined in (p7|). Using the values obtained in subsection 
III. A for the parameters uq and Eq entering , and choosing A such that the theoretical curves and the experimental 
data match reasonably, we find that z ~ 1 corresponds to a field B c ~ 0.025T. This is shown on Fig.|] where we note 
that the experimentally observed behaviour shows indeed a crossover at fields of the order of 0.05T. Therefore our 
theoretical crossover field B c is bit too small. However the two different regimes observed experimentally can be 
clearly interpreted in the framework of our theory: they describe the two situations where the diamagnetic response 
is essentially due to solely one type of vortex line elements (thermal for small B and field induced for large B). This 
is already gratifying for a first approach. For larger fields the term log(z) would yield a magnetization that increases 
more slowly with B than the measured data, which rather show a linear B-dependence. This is reproduced by the 
second part of ( p8| ) which, in our procedure, arises from the vortex correlations between different layers and becomes 
relevant for large values of B. 

These considerations are illustrated on FigjH where the data from Ref.6 are shown together with the theoretical 
curves given by Eq.(p8|). Although the quantitative agreement is less spectacular than in the case of the zero-field 
susceptibility, it still allows to extract the values of the "magnetically active volume fraction" A (from the small B 
region) and of the anisotropy parameter 7 (from the large B region). For the four temperatures of interest we find 
that the choice A ~ 0.15 is the most satisfying. This is quite reasonable for a layered compound such as underdoped 
YBCO where superconductivity occurs only in copper oxide planes which represent only a small fraction of the unit 
cell. Concerning the anisotropy 7, we find a value ~ 2. Its order of magnitude is correct since, by multiplying it by 
the lattice anisotropy d/a ~ 3, we obtain a value of 6 for the effective anisotropy. The latter is somewhat lower than 
what is observed in penetration depth measurements where it is found to be of the order of 25 at T c .tB However, we 
recall that 7 entered our theory through the vortex line effective stiffness in Eq.(|2(]). This rather rough description of 
the vortex loops interlayer correlations and the assumption we made by keeping independent of the temperature 
may possibly explain the fact that the observed value is a bit different from 25. 

Using the values of A and 7 discussed above, it is now possible to show that the second term in our expressions j2^ ) 
and ( |28| ) is irrelevant for the low-field behaviour, as we have anticipated when fitting the zero-field susceptibility. 
Finally we note that the picture of the vortex structure emerging from the above analysis which consists in loops and 
lines loosing their interplanar correlations over a very short distance n^d agrees very well with the usual descriptions 
of a n vor . t e x-system at temperatures lying above T c and above the vortex melting and a possible vortex decoupling 
lineflBEfB 

However, as Fig.^ shows, the experimentally observed crossover is extremely sharp - in fact much sharper than what 
we can obtain using our expression ((2^). This observation calls for a more refined treatment where other effects of 
the magnetic field than the mere creation of vortex lines must be taken into account. 

It is in fact possible to understand, in the framework of our approach, what makes the crossover so sharp. In terms 
of the susceptibility \(T,B) in Eq. (p3|) , the sharpness of the crossover suggests that the first term must tend to zero 
very rapidly as B increases. In this way, only the second (approximately constant) term remains and then yields the 
linear behaviour of the magnetization M(T, B) for B > B c . 

This behaviour could be achieved by assuming a small magnetic field dependence of the activation energy Eq and 
be motivated as follows. In Ref.21 it is emphasized that, below T c , the effective interaction between vortex loops 
is screened by the thermal defects of the Abrikosov vortex lattice. For small B, this effect is enhanced when the 
magnetic field is increased because the density of vortex lines is directly proportional to B as in Eq.(|27|). Above a 
crossover field B* the screening becomes weaker due to the fact that the finite stiffness (oc B 2 ) of the vortex lines 
inhibits the formation of further defects. In the liquid phase above T c which we are interested in, a similar qualitative 
behaviour can be expected, at least for low fields. A screened interaction between vortex loops then reduces the 
energetic cost of creating such an object. Thus the total number of thermally excited vortex loops must increase 
with the magnetic field B. Among them a (small) percentage corresponds to those contributing to the diamagnetic 
response as discussed in the previous section. They will also follow the above behaviour so that we may reasonably 
assume that their effective activation energy Eq introduced in Eq. (^5|) must decrease slightly when the magnetic field 
B increases. To lowest order, we have 




(29) 
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The large B behaviour is not important since it affects only the first term of the susceptibility when it is already 
very small. To illustrate this idea, we take a = 1 such that the value of the activation energy saturates at a value 
Ei 20 percents lower than E$ around B = 0.5T. After having performed a numerical integration of the susceptibility 
x(T, B) over B, we obtain the curves shown in Fig.^j. The crossover between the two field regimes is sharper than in 
the previous case and the corresponding value of the magnetic field B c is higher. 

We emphasize however that the above strategy was only intended to show how the developed formalism could be 
modified in order to improve the interpretation of the experimental results from Ref.6. A proper justification of the 
above assumption (^9|) requires concepts and methods which are beyond the scope of this work. 



IV. SUMMARY AND CONCLUSIONS 



We have derived explicit expressions for the field and temperature dependent diamagnetic susceptibility x(T, B) 
and magnetization M(T, B) of an anisotropic superconductor above its critical temperature. The superconducting 
fluctuations above T c are treated in the framework of the Lawrence-Doniach model including a magnetic field B 
perpendicular to the lattice planes. In order to describe specifically the precursor effects above T c , we have used the 
London approximation, assuming thereby that the relevant fluctuations are given by the phase of the order parameter 
only. Such an approach should apply to the underdoped regime of cuprates. In this context the current-current 
correlation function that determines x an d M is expressed by the structure factor of the relevant phase excitations 
- the thermally excited vortex loops and the field induced vortex lines. Our expressions for x an d M contain still 
several undetermined parameters such as the areal density of vortex line elements within a given layer, an appropriate 
length describing vortex line correlations between different layers, the anisotropy and a factor smaller than one that 
gives the fraction of the volume of the superconductor that is "active" in contributing to the fluctuation induced 
diamagnetism. These quantities, which must be compatible to another, are estimated by comparing x(T, B = 0) and 
M(T, B) to experimental data obtained on underdoped YBCOHand allow then to deduce an accurate picture of the 
phase system of such materials. 
Our findings can be summarized as follows: 

(i) There have been attempts to understand the same data based on Gaussian fluctuations of the anisotropic LD- 
model.H The corresponding fit for optimally doped YBCO, where the precursor region is much narrower, has worked 
out successfully. However, the discrepancies between experiment and the corresponding theoretical curves for the 
underdoped compound are very large, both in the temperature and magnetic field dependences. Our approach, based 
on topological phase excitations in a non critical regime, is much more satisfactory. 

(ii) The experimental low-field susceptibility has an activated temperature behaviour. This observation indeed points 
to excitations which have a finite creation energy (the vortex loops), rather than some wave- like fluctuations of the 
pairing field (the Gaussian modes). The value of the activation energy is roughly 20 times the critical temperature, 
what shows that the relevant vortex loops extend over several layers. The density of those loops is relatively small, 
indicating that, close to the critical temperature, the phase system consists predominantly of vortex loops running 
essentially parallel to the planes, and the line elements in z-direction which contribute to \ only develop significantly 
above the critical temperature. 

(Hi) The interlayer current correlations reflecting the geometrical structure of vortex loops in the z-direction are taken 
into account in terms of an effective correlation length £3 = n^d. The value of the latter has been determined by the 
vortex loop activation energy. In the temperature range of the measurements it is small, typically 2 or 3 interlayer 
distances. This feature can be also characterized by the stiffness of the vortex line which is found to be approximately 
inversely proportional to the anisotropy parameter. These observations indicate that the observed phase structure 
corresponds very well to the picture of weakly coupled layers of two-dimensional "pancakes" vortices. 

(iv) In the calculation of the molar susceptibility we have multiplied the molar volume by an overall prefactor 
called "effective active volume fraction" taking into account only those parts of the unit cell that contribute to 
superconducting fluctuations. This quantity was found to be of the order of 10-20 percents. 

(v) Both the activated behaviour of x and the small value of the above correlation length £3 — n^d point to the fact 
that the measurements of Ref.6 cover a "precritical" regime. Indeed, although the relevant physics is governed by 
phase fluctuations according to our assumptions, no divergence of the 3D XY-type has been observed. The latter is 
nevertheless expected for the data taken the closest to T c , as it has been identified in Ref.8 for layered superconductors 
with a reasonable value of the anisotropy. In fact one observes that some deviations from the activated behaviour 
appear for the temperatures the closest to T c . This may indicate a crossover to a critical regime in which the In- 
dependence of x is n0 more simply governed by the total vortex density, but by a diverging correlation length £ 
associated with the free vortex density which should be used instead of nv in the structure factor ( |l4j ) and in x 
[Eq.(23)], as it was done in the 2D case for calculating the zero field susceptibilityllj Moreover, in this case the 
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dimensionless quantity z defined in Eq.(p7|) becomes the scaling variable -B£ 2 /$o which has been used extensively to 
study critical properties under magnetic fieldJlj 

(vi) The experimental data show an extremely sharp crossover in the field dependence of the magnetization. Within 
our approach this is attributed to a subtle interplay between the two types of vortices suggesting that the presence of 
field-induced vortex lines is favoring the creation of thermal vortex loops containing at least two segments along the 
^-direction. However, even our most refined theoretical curves are still more smooth than the data. It is therefore not 
excluded that the observed behaviour may have yet another origin than superconducting fluctuations, for instanced 
sample inhomogeneities, as it has been suggested and observed in other samples by the authors of the experimental 
workQ that we have analysed. 

We thank X. Zotos, A. A. Varlamov, S.R. Shenoy and T. Schneider for interesting discussions. This work has been 
supported by the Swiss National Science Foundation (projects 20-49586.96 and 2000-053697.98/1). 
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APPENDIX: 



Here we show that the zero wave vector value of the positional structure factor ([12|) compensates the singular term 
of the current-current correlation function ([ll]) . 
Starting from Eq.(|l2|), we find 

5(q = 0)= t(s,n)t(s',n') 

= J2( N +- N -) 2 

where N + [N-] is the number of vortex line elements with t(s,n) = +1 [t(s,n) = —1] in layer n. They correspond 
exactly to the 2D concepts of vortex and antivortex which are often interpreted as positive and negative charges of a 
2D Coulomb gas. 

The external flux generates Np "charges" according to 

N F = (A2) 

Due to the incompressibility of the magnetic field (V-B = 0), this number is the same for all layers. 

The neutrality condition requires that, in every layer n, the external flux is exactly compensated by the net charge 

iV_|_ — of the vortex-antivortex system: 

N F = N+ - N_, Vn. (A3) 



Therefore (Al) may be written as 



^(q = 0)=A 2 (^) 2 (A4) 



which is exactly the singular term of the current-current correlation function C (q) . 
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FIG. 1. Schematic representation of a 3D thermally excited vortex loop (left) and of a field-induced vortex line (right). Since 
in the Lawrence-Doniach approach the phase field is defined only for discrete values of the ^-coordinate (z = nd), one often 
refers to these structures as "stacks of pancake vortices". The latter are represented by the grey ellipses whereas the lines 
linking them are then just guides to the eye. 
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FIG. 2. Zero-field susceptibility \(T, B = 0) : the triangles (a) are the experimental data from Ref.6 and the full line is the 
best fit of the theoretical expression (B3J) assuming the activated behaviour (123). 
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FIG. 3. Magnetization M(T — const, B) : the full lines correspond to the theoretical expression ( |28| ) with the values of the 
parameters mentioned in the text and the triangles (a) in the inset are the experimental data from Ref.6. 
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FIG. 4. Improved version of the magnetization M(T — const, B) : the full lines correspond to the theoretical expression 
based on Eq.([29|). 
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